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INTRODUCTION 


I  . 

A .  BACKGROUND 

Historically,  the  attitude  stability  conditions  of 
spinning  spacecraft  have  been  derived  using  the  energy  sink 
method.  A  significant  shortcoming  of  this  method  is  that  it 
does  not  account  for  the  dynamic  interaction  of  liquid  motion. 
Tor  current  spacecraft  with  large  amounts  of  liquid 
propellant,  dynamic  interaction  can  destabilize  the 
spacecraft.  As  a  result,  under  some  circumstances  the 
stability  prediction  of  the  energy  sink  method  can  be  in 
error . 

Despite  this  shortcoming,  the  energy  sink  method  continues 
tc  be  in  wide  spread  use  as  an  analytical  technique  fcr 
determining  attitude  stability  for  dual  spin  spacecraft.  This 
being  the  case,  it  is  essential  thi1  the  nature  of  the 
deficiency  and  the  conditions  under  which  it  occurs  be 
understood . 

B.  OBJECTIVE 

The  objective  of  this  study  is  to  investigate  and  more 
accurately  define  the  nature  of  this  shortcoming. 
Specifically,  the  case  of  a  marginally  stable,  dual  spin 
spacecraft  with  a  despun  platform  will  be  explored.  The 


stabili1 


the  spacecraft  will  be  examined  a: 


•  a  ■■■  o  v*  ^  v 


.  c  ^  «  4-  r  ^  , 


asymmetry  ar.d  the  fuel  load  are  varied. 


Marginally  stable,  in  terms  of  the  energy  sink  stability 
conditions,  means  the  spacecraft  has  an  inertia  ratio  slightly 
greater  than  one.  The  percent  asymmetry  is  defined  as  the 
ratio  of  the  difference  of  the  two  transverse  moments  of 
inertia  to  the  sum  of  the  two  transverse  moments  of  inertia 
times  100.  The  fuel  load  is  the  amount  of  fuel  on  beard 
expressed  as  a  fraction  of  the  total  fuel  capacity  cf  the 
spacecraft . 

To  perform,  the  analysis,  a  computer  simulation  developed 
by  Chung  [Ref.  1]  of  a  dual  spin  spacecraft  will  be  used  to 
determine  the  stability  of  the  various  conf igurati ons .  These 
results  will  be  compared  to  the  predictions  cf  the  energy  sink 
stability  criteria  as  developed  by  Likins  in  Reference  2. 


C.  LITERATURE  REVIEW 

The  energy  sink  approach  was  first  applied  to  study  the 
effects  of  energy  dissipation  on  the  stability  of  a  freely 
spinning  body  in  1963  by  Thomson  and  Reiter  [Ref.  3],  This 
led  to  the  well  known  requirement  for  stability,  an  object 
must  spin  about  its  axis  with  the  largest  principal  moment  of 
inertia . 

In  1966,  Likins  [Ref.  2]  developed  the  energy  sink 
stability  conditions  for  asymmetric,  dual  spin  spacecraft 
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to  the  average  of  the  rate  of  change  of  the  corresponding 
inertial  so  in  sate.  He  sen t  on  to  ocint  cat  list  the  ratio 
of  the  spin  moment  of  inertia  to  the  algebraic  mean  transverse 
moment  of  inertia  was  the  critical  stability  quantity. 

In  1972,  two  papers  questioned  Likins'  conclusions 

"  V  w  ”  ''  3  ^  "J  1*-  3 't.  **  w  0  mj  ^  *•  * 7  tl  31  "  ^  o  ■  ;  c  -  *-  V-  ^  *•  -  —  ->  ■*- 


platform. 

In  1983,  Cochran  and  Shu  [Ref.  6]  used  the  generalized 
method  of  averaging  to  study  both  energy  dissipation  and  the 
energy  addition  required  to  maintain  the  constant  spin  rate 
of  the  rotor.  Their  conclusions  substantiated  Hubert's 
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hypothesis  concerning  rore  energy  assuming  that  the  internal 
mass  motion  is  sufficiently  small. 

In  Reference  1,  Chung  studied  the  application  of  the 
energy  sink  method  to  the  INTELSAT  VI  satellite.  His 
conclusion  was  that  the  energy  sink  method  did  not  correctly 
predict  stability  for  all  cases.  The  results  of  Reference  7, 
also  a  study  of  the  INTELSAT  VI  satellite,  indicated  that 
stability  increased  as  the  fuel  load  increased,  and  that 
stability  decreased  as  the  platform  asymmetry  increased. 

The  energy  sink  method  continues  to  be  used  throughout  the 
community  to  determine  the  stability  of  spinning  spacecraft. 
In  order  to  have  confidence  in  the  predictions  of  the  method, 
if  is  necessary  to  understand  its  limitations.  It  is  hoped 
that  this  study  will  provide  additional  insight  into  the 
energy  sink  method's  transition  cone.  In  other  words,  where 
the  boundary  is  between  accurate  energy  sink  predictions  and 
erroneous  predictions. 

D.  ORGANIZATION  OF  STUDY 

Chapter  II  presents  the  derivation  of  the  energy  sink 
attitude  stability  criteria  for  an  asymmetric  dual  spin 
spacecraft,  a  description  of  the  simulation  used  to  model  the 
satellite  motion,  and  a  description  of  the  satellite 
configuration  used  for  the  study.  Chapter  III  describes  the 
development  of  the  system  parameters  for  the  spacecraft 
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With  the  coordinate  system  located  at  the  rents:-  of  "".ass  cf 
the  total  system,  and  fixed  in  the  body  P  such  that  the  center 
cf  mass  cf  both  bodies  lies  along  the  03  axis,  the  rotational 
equation  cf  motion  governing  the  system  is, 

M  =  dH/dt  =0  '1) 
where  M  is  the  total  external  moment  exerted  on  the  system  and 
H  is  the  total  angular  momentum  of  the  system. 

The  resulting  linearised  equations  c f  motion  are. 


M ■  -  Z  ^ 

'  T  ,  -  I  '  '  w  -  W  ' 

-  r  '(7  -  1 

'  ->  \ 

M-  =  I-w-  J- 

(T:  -  I.lw.’j' 

-  :rw..w-  = 

a 

f  31 

M-  =  T.w.  * 

I  w„  *■  I-w  w- 

■”  ^  •  V  •  'n  *  — 

0 

'  4) 

where  I.  is  the  total  pri 

n c i p  a  1  axial 

moment  of 

inertia 

of  the 

2  "3 1  ,  7  is  principal  axial  rr. c m s n f  cf  i  r. a r f. i a  cf  flc© 

body  ?. ,  I-  and  I-  are  the  principal  transverse  moments  cf 
inertia  of  ^be  s v s ^  is  *■  *■“ 0  ax'  al  * c nr  c n e n t  cf  tbe 

angular  velocity  of  the  body  P,  w,  is  the  relative  rate  at 
which  the  body  R  is  moving  with  respect  to  the  body  ?,  and  w-_ 
and  w-  are  the  transverse  components  of  the  total  system’s 
angular  velocity.  These  three  equations  have  four  unknowns: 
vr,  w* ,  w.,  and  wr .  The  fourth  equation  describes  the 
rotor/platform  interface, 

=  Ir(wp  +  wr)  (5) 

where  T„  is  the  resultant  moment  of  all  the  forces  acting  by 
P  on  R  about  the  rotor  axis.  Equations  (2)  and  (3)  can  be 
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0 


■  =5  w  r  i 


as  , 


W-  *■  3-W»  = 


W"  -  0-w-  =  0 


where , 


(  ) 
(7) 


a ■  =  [  ( I. 

-  :-)w. 

4- 

I..Wr]/I- 

'  3 ) 

o;  -  [(Ir 

-  ::)w. 

+ 

Irwr]/x: 

(9) 

Letting  w-(0)  =  w*  and  w-(0) 

•H 

o 

II 

e  1 

ds  the  following 

solution 

t o  equations  (61  and  ( ”  )  , 

w-  =  -  w-v"  ( c  •  /a- )  sin."  (  c  ■  o- '  t  (10) 


'.'p'.yi-yq  the  ?.cut  h-Hurw  1 1  z  stability  oriteri 
characteristic  equation  of  equations  ( •' )  and  ( 7  N 


l'*.  *  T  V.  ~  ±» 

.  .  U  ^  Z)  -  . .  c*  - 


for  s t abi 1 i t y , 

o-o*  >  0  ( 12  ) 

The  equations  for  the  angular  momentum  and  rotational 
kinetic  energy  for  the  system  described  above  are, 

h-  =  I-~wp-  +  I-’w-"  +(  I.w.  +  Irwr) ■  (13) 

2T  =  I-W--  +  I^wi*  +  Ir-w--  +  Irw„‘  +  2I:wrw-  (14) 

The  nominal  angular  momentum  is  given  by, 

hr  =  IpWp  +  Irwr  (15) 

and  is  assumed  to  be  constant.  Taking  the  derivative  of 
equations  (13)  and  (14)  with  respect  to  time,  and  taking  into 
consideration  conservation  of  angular  momentum  and  assuming 
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an  energy  dissipation  mechanism  yields, 
2  =  I-  w-w-  I-'W-W-  x  (I.w.  x  I  ..W 


m  =  l-w-w-  I-W-W-v  +  I„w.w„  +  Lw/<j,  + 


+  I^wrwr  +  <  0  (17) 
New  the  solutions  in  equations.  (10)  and  (11)  are  no  longer 
correct.  However,  if  the  effects  of  the  energy  dissipation 


mechanism  are  felt  slowly,  a  solution  of  the  same  f 


"ay 


assumed  but  with  w*  =  w-.(t),  a  slowly  varying  function  of  time. 


Substituting  this  new  solution  into  equations  (Id)  and  'l7), 


the  period  t  -  2 r./ o*,r(  o- / c- ) 


»  1  i  m  '  n  o  ^  A  q  pj  ^  ^ 


w«w«  =  -2o-h->(  Icw:  x  I^Wv)/(I-'C'  *■  I-'o-) 

T  =  W'W"(I’0'  +  I-o*  )/2o-  +  I-w.w,.  x  I  ,y,w,.  + 


+  I„w,.wc  +  irw,w. 


(19) 


Substituting  equation  (18)  into  equation  (19)  and  introducing, 
o*  =  h:(I:o:  +  I:o-)/(I::o-  +  I;-o:)  (29) 

with  a  little  manipulation  produces. 


f  =  -(I-  -  Ir)ws(c-  -  wt)  x 

-  Ir(wp  +  Wr)(o;  -  (w3  x  wr)) 


(21) 


=  Pp  +  pr 


where  P?  and  Pr  are  the  platform  and  rotor  components  of  the 
energy  dissipation  rate.  It  is  necessary  that, 

PP  <  0 


P,  <  0 
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:inw.  let 


r;  -  O'  -  w. 

t  ->  o  \ 

Or  =  o*  -  (w-  +  w_) 

(23) 

so  that, 

?-  =  -( Is  - 

(24) 

?.  =  -I,(w.  +  Wr)or 

(25) 

or  alternately 

p./o-  =  -(Ir  -  I:)w. 

(26) 

?r/o..  =  -  wr) 

(  O'7  ) 

Substituting  equations  (26)  and  ( 2? )  into  equation 

'  3 ) 

yields , 

W  ~  W  p  ~ 

[2o:h:/(I::o:  t  I:'o:)  ]  [  (P./Cp)  -  (?r/or)] 

(23) 

As  a  necessary 

and  sufficient  condition  for  stability, 

w-w-  <  0 

(29) 

Since  >  ^ 

is  also  necessary  for  stability  and  hq  > 

C  by 

convention,  it 

follows  that. 

(P./O.)  A  (Pr/Or)  <  0 

F  t  “ 

(30) 

And  because  P.  and  P„  are  both  negative,  at  least  one  of  the 
two  of  o.  and  o,  must  be  positive  such  that  the  total  effective 
energy  dissipation  rate  of  equation  (30)  is  negative. 


In  terms  of  the  specific  application  to  this  study  the 
equations  simplify  as  follows.  For  a  system  with  a  despun 
platform  and  no  damping  mechanism  on  the  platform,  w?  is 
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essentially  zero  and  P„  is  zero.  Therefore, 

°1  =  1  :w:// 1 : 
o,  =  Irwr/I, 

and , 

O;  =  IrMr(Il°:  +  I:o,)/(I,:o1  +  I.-o:) 
hn  =  I.w. 

Now , 


'3) 

(9) 

(20) 


(15) 


o.  =  o.  (22) 

3,  =  0;  "  Wr  (23) 

and  the  stability  criteria  is, 

pr/o.  <  o 

Or,  since  P„  must  be  negative, 

or  >  0  '?!) 

Substituting  equations  (8)  and  (9)  into  equation  (20)  yields, 

o„  =  2I„w,/  ( I,  +  I,)  (32) 

*.  *.  * 

and  substituting  equation  (32)  into  equation  (23)  gives, 

or  =  2I:wr/(I,  +  I,)  -  wr  (33) 

or,  from  equation  (31), 

wr[2Ir/(I,  +  I,)  -  1]  >  0  (34) 

Dividing  by  w„  and  adding  one  to  both  sides  of  equation  (34) 

yields , 


217(1,  +  I2)  >  1  (35) 

which  is  the  familiar  spin  axis  moment  of  inertia  divided  by 
the  algebraic  mean  transverse  moment  of  inertia  developed  by 
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Likins.  It  is  at  this  point  that  Spencer  asserts  that  using 
the  geometric  mean  in  the  denominator  yields  more  accurate 


r esul ts . 


Spencer's  conclusion  is  based  on  the  fact  that 
rotor-fixed  nutation  frequency,  actually  varies  over 
a  function  of  I.  ,  the  transverse  moment  of  inertia, 
states  that  I-.  varies  as, 

T.(wt)  =  I  •  sin'wt  +•  I-ccs'wt 


w ’  ,  the 
time  as 
And  he 

(36) 


where . 


w  -  :  ( :•:«,)  ]w„ 

quaticns  'IS'  and  (11).  These  variations 


(37) 

'  m  V-, 


reflected  in  T.  Likins  neglected  this  in  averaging  over  the 
period  to  obtain  equations  (IS)  and  (19).  Unf  or  t  ur.at  el  y ,  this 
makes  the  solution  for  an  asymmetric  satellite  much  more 
complicated.  Spencer  does  not  provide  a  complete  solution. 
Instead,  he  presents  a  simplified  example  and  shows  that  using 
the  geometric  mean  for  I*-  provides  a  closer  approximation  to 
simulation  results. 


B.  SIMULATION  DESCRIPTION 

The  computer  simulation  used  in  this  study  to  determine 
the  numerical  solution  to  the  equations  of  motion  of  the 
satellite  was  developed  as  a  part  of  the  work  done  by  Chung 
in  Reference  1.  This  section  provides  a  description  of  this 
simulation. 
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A  dual  spin  spacecraft  can  be  modelled  as  two  rigid  bodies 
capable  of  rotating  relative  to  each  other  about  a  common 
axis.  This  common  axis  passes  through  the  mass  center  of  the 
body  representing  the  rotor.  The  liquid  in  spherical  fuel 
tanks  can  be  modelled  as  an  axisymmetric  spherical  pendulum 
with  the  hinge  point  at  the  center  of  the  tank.  The  spherical 
pendula  can  be  mounted  on  either  of  the  two  bodies  that 
constitute  the  system  (Figure  2).  The  energy  dissipation  due 
to  the  liquid  sloshing  can  be  included  in  the  model  by 
introducing  viscous  damping  in  the  spherical  joint  of  the 
pendul a . 


Figure  2.  Schematic  of  a  generic  dual  spin  spacecraft. 
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Chung  applied  Kane's  method  to  the  model  described  above 
to  determine  the  equations  governing  the  dynamic  behavior  of 
the  system.  This  method  entails  defining  generalized  forces 
called  inertia  forces  and  active  forces.  These  forces  are 
expressed  as  functions  of  the  generalized  speeds  (rotational 
velocities)  of  the  various  components  of  the  system.  The 
total  force  acting  on  the  system  can  be  summarized  as, 


=  0  (r  =  1,  2,  3 . N) 


'8 ) 


where  N  is  the  total  number  of  generalized  speeds,  and  and 


^  P  v  o  ** 


the  total  generalized  active  and  inertia  forces  in 


i.  n  0  it  Z-cil  sr>ac e  ^^0  ^ v  /-> ® ^y 

is  equal  to  the  number  of  degrees  of  freedom  of  the  system. 
In  this  case,  three  for  the  three  components  of  the  angular 
velocity  of  the  platform,  plus  one  for  the  relative  velocity 
between  the  platform  and  the  rotor,  plus  three  times  the 
number  of  fuel  tanks  (each  pendulum  representing  a  fuel  tank 
has  three  degrees  of  freedom) . 

To  obtain  the  generalized  inertia  forces  of  the  system, 
the  contributions  from  the  platform,  rotor,  and  pendula  are 
summed.  To  obtain  the  generalized  active  forces  of  the 
system,  the  contributions  from  all  active  forces  on  each  part 
of  the  system  are  summed.  It  is  assumed  that  the  resultant 
of  the  external  forces  of  the  platform  and  the  rotor  are  zero 
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and  that  no  external  forces  act  on  the  pendula  mounted  on 
either  the  platform  or  the  rotor. 

In  addition  to  the  generalised  speeds  obtained  by  solving 
the  differential  equations  developed  using  Kane's  method, 
several  other  quantities  are  useful  in  understanding  the 
motion  of  the  spacecraft.  These  include  the  central  angular 
momentum  of  the  system,  the  kinetic  energy  of  the  system,  the 
energy  dissipated  through  the  spherical  joints  of  the  pendula, 
the  work  done  by  the  motor  and  the  external  forces,  and  the 
nutation  angle  of  the  system.  All  of  these  can  be  expressed 
in  terms  of  the  generalized  speeds. 

The  simulation  takes  as  input  the  system  parameters  that 
characterize  the  properties  of  the  satellite  being  simulated. 
These  include  the  mass,  moments  of  inertia,  location  and 
orientation  for  each  component  of  the  system  relative  to  the 
center  of  mass,  as  well  as  the  key  properties  of  the  pendula 
(length  and  damping  coefficient).  The  initial  conditions  for 
all  of  the  generalized  speeds  and  coordinates  must  also  be 
provided.  The  output  of  the  simulation  is  a  set  of  values 
that  characterize  the  state  of  the  system  at  a  given  point  in 
time.  These  include  the  nutation  angle,  the  kinetic  energy, 
the  total  energy,  the  work  done  by  the  rotor  and  the  external 
forces,  the  two  components  of  the  transverse  angular  velocity, 
the  platform  angular  velocity  and  the  rotor  angular  velocity. 
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C.  SATELLITE  CONFIGURATION 

The  configuration  of  the  satellite  used  in  this  study  is 
a  derivative  of  the  INTELSAT  VI  satellite.  The  satellite 
consists  of  a  platform  and  a  symmetric  rotor.  There  are  four 
fuel  tanks  and  four  oxidizer  tanks  mounted  on  the  platform  in 
the  configuration  shown  in  Figure  3. 


Figure  3.  Arrangement  of  fuel  and  oxidizer  tanks. 
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These  parameters  are  modified  using  the  procedure  described 
in  Chapter  III,  Section  A  to  obtain  the  configuration  used  in 
the  study  (Tables  II-IV). 
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PROCEDURE 


This  chapter  describes  the  methodology  used  ever  the 
course  of  the  study.  The  first  section  deals  with  how  to 
determine  the  spacecraft  ccr.f i  ~urat  i  on  and  system  parameters. 
Section  3  describes  how  the  energy  sink  criteria  is  applied 

■■ "/ r.  f  i.  ^ "U r  3  -  1.  o n  3  c  ^  i  o n  C  c  o ^  ti  3  i.  s  u  ^  ^  r  cr  c  —  •  •  -  ^  ^  ~ 

the  simulation  ar.d  obtaining  the  simulation  output. 

A.  SYSTEM  PARAMETERS 

Tr.  order  to  investigate  the  behavior  of  a  marginally 
stable  spacecraft  it  is  first  necessary  to  model  a  spacecraft 
with  this  configuration.  To  accomplish  this,  .he  system 
parameters  of  the  INTELSAT  VI  satellite  were  modified  to 
obtain  the  appropriate  configuration.  Cne  of  the  programs 
written  by  Chung  in  Reference  1  to  support  the  simulation 
produces  as  its  output  a  summary  of  the  INTELSAT  VI  system 
parameters.  To  modify  the  system  parameters  to  achieve  a 
given  inertia  ratio,  l./h,  where  I  =  (!'  +  I-)/2,  the 

algebraic  mean  of  the  spacecraft  transverse  moments  of  inertia 
must  be  calculated.  This  is  then  multiplied  by  the  inertia 
ratio  to  determine  the  desired  axial  moment  of  inertia,  Is. 
Per  the  case  of  a  despur.  platform,  this  equates  to  the  axial 
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moment  c  f  :r.  artia  of  the  wet  rotor,  I..  However,  the 
siT.ulaticr-  takes  as  its  input  the  dry  spacecraft  parameters. 
Subtracting  the  old  wet  rotor  moment  of  inertia  from  the 
desired  moment  of  inertia  yields  the  amount  the  dry  rotor 
axial  moment  of  inertia  must  be  increased.  It  is  a  general 
theorem  of  rigid  body  mechanics  that  for  a  given  body,  the  sum 
:f  ar.y  two  ;f  the  principle  moments  of  inertia  must  be  greater 
than  the  third.  This  must  be  kept  in  mind  to  achieve  a 
realistic  design .  To  adhere  to  this  principle,  the  sum  of  the 


dry  rite:  transverse  moments  ct  inertia  is  s: 
new  dry  r'tcr  axial  moment  of  inertia.  The  result  is  divided 
by  two  and  added  to  each  of  the  dry  rotor  transverse  moments 
of  inertia.  An  addi t i cnal  adjustment  is  made  to  achieve 
perfect  symmetry  on  the  rotor.  The  difference  between  the  dry 
rotor  transverse  moments  of  inertia  is  added  to  the  smaller 
of  the  two  to  make  them,  equal.  To  balance  the  system,  and 
maintain  the  same  total  spacecraft  transverse  moments  of 
inertia,  the  corresponding  platform  transverse  moments  of 
inertia  are  decremented  by  the  same  amounts  as  were  added  to 
the  rotor  transverse  moments  of  inertia. 


B.  ENERGY  SINK  PREDICTION 

Given  numerical  values  for  1^,  I;-  Ip/  Ir,  w5,  and  wr, 
developing  an  energy  sink  prediction  is  a  reasonably  straight 
forward  application  of  the  equations  derived  in  Chapter  II. 
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and  that  P.  and  ?,  are  always  negative.  Concentrating  on  the 
second  relationship,  for  the  system  being  studied,  ?.  is  zero 
and  ?.,  is  negative  since  all  of  the  fuel  tanks  are  mounted  or. 
the  rotor.  Therefore,  if  o,  is  positive  the  system  is 
inherently  stable.  However,  if  o,  is  negative  o„  must  be 
positive  and  a  nutation  damper  must  be  installed  on  the 
platform  to  overcome  the  destabilizing  effect  of  the  rotor. 

C.  SIMULATION  PROCEDURE 

To  determine  the  stability  of  the  system  the  parameter  of 
interest  is  the  nutation  angle,  5.  If  $  grows  without  bound, 
the  system  is  unstable.  If  $  damps  out  to  zero,  the  system 
is  stable.  In  executing  the  simulation,  an  initial  transverse 
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nutation  ingle. 

The  first  step  in  the  study  entails  finding  a  satellite 


configuration  with  an  inertia  ratio  slig 


^  ..UJ-  1  ~  *“  V. 


eater  tr.ar.  one 


u  v  •*  v  *■ v 


.aticn  results  indicate  the  svstem  is 


stable.  An  inertia  ratio  of  1.C1  is  arbitrarily  chosen  for 
the  initial  run.  If  the  simulation  results  indicate  the 
system,  is  unstable  the  inertia  ratio  must  be  increased  and  the 


simulation  run  again.  This  orocess  is  continue 


w  .-t-i  1 


table  configuration  is  found  up 


n  1  a  -T”  ^  '  r* ' '  ^  -n  *- 


z  t  *  \  d  v  will  be  based.  After  let  e  r  tt*.  i  n  i  n  cf  3.  stable  sonfi  juration 
the  platform  asymmetry  is  varied  and  the  satellite  metier,  is 
simulated  for  three  different  fuel  loads.  The  INTELSAT  VI 
satellite  uses  a  liquid  propellant  rocket  for  its  apogee  kick 
motor  and  ccnsum.es  nearly  75%  of  its  fuel  during  this 
maneuver.  The  beginning  of  life  fuel  load  for  the  INTELSAT 
VI  satellite  upon  which  the  study  configuration  is  based  is 
25.2%.  Therefore,  fuel  loads  of  26.2%,  20%.  and  15%  were 
chosen  for  this  study.  The  initial  asymmetry  was  arbitrarily 
selected  to  be  5%.  If  the  results  indicate  the  system  is 
stable  the  asymmetry  is  increased  by  5%  and  the  simulation  run 
again.  This  process  continues  until  the  simulation  results 
indicate  the  system  is  unstable. 


-a  uS*  -  ^  w  : 
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RESULTS  AND  ANALYSIS 


This  chapter  presents  the  results  of  the  study  and  a 
discussion  of  their  s igni f i cance .  The  first  section  presents 
a  summary  of  the  results  of  the  simulation  runs.  Section  5 
contains  a  summary  of  the  energy  sink  predictions  for  all  of 
the  cases  simulated. 

A.  SIMULATION  RESULTS 

As  outlined  in  chapter  ITT.  the  first  steo  in  the  study 
was  to  find  a  configuration  for  which  the  system  was  stable. 
The  first  simulation  run  was  for  the  case  of  a  symmetric 
system  with  an  inertia  ratio  cf  1.01  and  a  fuel  load  of  26.2%. 
The  simulation  time  was  set  at  200  seconds.  The  result  is 
shown  in  Figure  4. 


on  cnB:T,  SYM,  !*>/  >t  ■-  i  oc,  ?s  cufl 


TiMECXQ 

Figure  4.  Initial  simulation  result. 
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From  the  graph  it  is  difficult  to  determine  whether  the  system 
is  stable  or  unstable  because  the  simulation  time  was  not  of 
sufficient  length,  r.ather  than  extending  the  s  imu!  a*-  i  nn  t  i  me 
to  arrive  at  a  more  conclusive  result,  the  simulation  was  run 
again  with  the  damping  coefficients  increased  by  a  factor  of 
100.  Figure  5  shows  that  this  time  the  results  are  much  more 
definitive . 

UN  ORB  I  T  ,  bYM  ,  I  S/  It  —  1  .  G  I ,  <^u  .  I’  uhi_ 


JUU*’  -1UU1 


Figure  5.  Initial  simulation  result  with  increased 
damping . 

The  system  is  obviously  unstable.  From  this  point  on,  all 
simulation  runs  were  made  with  the  damping  coefficients 
increased  in  order  to  minimize  the  main  frame  CPU  time. 
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Next  the  simulation  was  run  with  an  inertia  ratio  of  1.1. 
The  result  reveals  that  the  system,  is  stable  as  shown  in 
Figure  6 . 

ON  ORB  I T  ,  SYM,  is/  It  =  1.1,  26  2%  FUEL 


90022001 


0  10  20  JO  «0  50 

r  ifctcsec} 

Figure  6.  Result  for  an  inertia  ratio  of  1.1  and  a  fuel 
load  of  26.2%. 

Having  achieved  stability,  the  fuel  load  was  changed  to 
20%  and  the  simulation  was  run  again.  The  system  remained 
stable  as  shown  in  Figure  7  so  the  fuel  load  was  reduced  to 
15%  and  the  simulation  run  again.  The  result  still  indicated 
that  the  system  was  stable  (Figure  8).  Having  achieved 
stability  for  all  three  fuel  loads,  the  next  step  was  to  begin 
varying  the  asymmetry  of  the  platform. 
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90022003 
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Figure  7.  Result  for  an  inertia  ratio  of  1.1  ar.c 
load  of  20%. 


ON  ORB i T ,  SYM,  I s/ i t  =  1.1,  15.0%  FoEL 


90022002 


Figure  3.  Result  for  an  inertia  ratio  of 
load  of  15%. 


1 . 1  and 
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A  summary  of  the  results  for  si’  the  ssy-^^etry 

variations  can  be  seen  in  Table  V.  The  system  remained  stable 
for  every  case  simulated. 

TABLE  V.  Simulation  results,  2 1  s  /  ( I  •>  +  Ii)  = 

1.1. 


FUEL  ASYMMETRY 

LOAD  5%  10%  15%  20%  25%  35%  55% 

15%  S  S  S  S  S  S  S 

20%  S  3  S  S  S  S  S 

26.2%  S  S  S  S  S  S  S 

Note:  U  =  Unstable,  S  =  Stable 


In  light  of  the  results  reported  in  Reference  7,  the 
results  in  Table  V  were  surprising.  However,  upon  re¬ 
examining  the  energy  sink  equations,  it  was  discovered  this 
is  exactly  what  should  have  been  expected.  Starting  with  a 
svmm.etric  satellite,  such  that  Ii  =  I»  =  I*,  the  stability 

"  A  —  - 

criteria  is , 

2Ir/(I:  +  I2)  >  1  (35) 

To  create  an  asymmetry,  let  Ii  =  It  +  a  and  I2  =  I*  -  a. 
Substituting  for  I2  and  I2  in  equation  (35)  reduces  to  the 
identical  equation  indicating  the  system  should  remain  stable. 
Using  the  geometric  mean  as  suggested  by  Spencer,  equation 
(35)  becomes, 
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▼  f  ->  O.  \ 

Substituting  fcr  I-  and  I-  here  yields, 

-  a' )  )  >  1  f  40  > 

which  has  to  be  greater  than  the  original  inertia  ratio 
indicating  that  the  system  is  getting  more  stable  as  the 
asymmetry  is  increased.  Regardless  of  the  method  used  to 
determine  the  average  transverse  moment  of  inertia,  both 
T?  ^  0  cl  1.  c  v  ♦■no  c;  vgt  S  h  C  U  1  ci  ^  1  ^  3  t  Hb  '  b6r5'aS3  1.  v-  /  1  iw 

g^oater  their*  c  n  0 .  This  b  0  i  n  9  th.0  0030.  it  was  new  necessary 
to  find  a  configuration  which  was  unstable  to  determine  if  the 
system  would  become  stable  as  the  asymmetry  was  increased. 

The  simulation  was  run  for  symmetric  configurations  with 
inertia  ratios  of  1.03,  1.05,  and  1.07.  The  results  are 

summarised  in  Table  VI.  They  shew  that  fcr  a  symmetric 
spacecraft  the  stability  cutoff  falls  between  inertia  ratios 
of  1.03  and  1.05. 

TABLE  VI.  Simulation  results  of  varying  inertia  ratio. 


FUEL 

LOAD 

1.01 

INERTIA  RAT 
1.03  1.05 

10 

1.07 

1 . 1 

15% 

- 

U 

S 

S 

s 

20% 

- 

U 

S 

s 

s 

26.2% 

u 

U 

s 

s 

s 

Note:  U  =  Unstable,  S  =  Stable 


(?■*  yom  ^  V  ^  ^  ra  ^  ^  o  ^/T  ■»«*>*■»  *^  ■»  ^*  V*  3r* 

inertia  ratio  of  1.03  was  selected  to  evaluate  the  influence 
of  varying  the  asymmetry. 

The  summary  of  these  results  is  provided  in  Table  VII.  As 
expected,  the  system  eventually  became  stable  as  the  asymmetry 
was  increased. 

TABLE  VII.  Simulation  results  for  varying  the 
platform  asymmetry. 


FUEL  ASYMMETRY 

LOAD  0%  25%  40%  55% 

15%  tj  U  S  S 

Q  ^  ry  r*  C  C* 

2  6.2%  U  U  S  S 

Mote:  U  =  Unstable.  S  =  Stable 


what  appeared  to  have  happened  was  that  as  the  asymmetry 
was  increased,  the  inertia  ratio  increased  enough  to  become 
greater  than  the  stability  cutoff  inertia  ratio  as  equation 
(40)  indicates.  To  confirm  this,  a  graph  of  inertia  ratio 
versus  asymmetry  was  created  using  equation  (39)  to  employ  the 
geometric  mean  transverse  moment  of  inertia  (Figure  9).  As 
expected,  it  showed  that  for  a  symmetric  inertia  ratio  of 
1.03,  as  the  platform  asymmetry  was  increased,  the  inertia 
ratio  increased  enough  to  become  greater  than  the  stability 
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cutoff  inertia  ratio  that  falls  between  1.03  and  1.05.  This 


plot  combined  with  Table  VII  indicate  that  the  stability- 
cutoff  inertia  ratio  falls  between  1.032  and  1.035. 
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Figure  9.  Inertia  ratio  vs.  asymmetry  for  26.2%,  20%, 
and  15%  fuel  loads. 


To  try  and  pin  down  the  stability  cutoff,  the  simulation 
was  run  for  the  symmetric  case,  for  all  three  fuel  loads  with 
inertia  ratios  of  1.035,  1.04  and  1.045.  Adding  the  results 
to  Table  VI  shows  that  the  stability  cutoff  inertia  ratio  for 
a  symmetric  satellite  is  different  than  for  an  asymmetric 
satellite.  The  simulation  results  indicate  the  stability 
cutoff  for  the  symmetric  satellite  is  approximately  1.045 
(Table  VIII) . 
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TABLE  VIII.  Simul  aticn  results  for  virvir.g  inertia  ratio 
of  symmetric  spacecraft. 


FUEL 

roin  ■» 

n  y 

1  a  ? 

INERTIA 
1.035  1.04 

RATIO 

1  A  ^  C  1 

A  C 

1  A  "7 

1  1  0 

15% 

- 

u 

U 

u 

M 

S 

C 

s 

20% 

- 

u 

u 

u 

M 

S 

s 

s 

25.2% 

u 

u 

u 

u 

M 

s 

C 

s 

Note:  U  =  Unstable,  S  =  Stable,  M  =  Marginal 


Finally,  several  simulation  runs  were  made  with  fuel  loads 
of  50%  and  75%  to  investigate  the  effect  of  higher  fuel  loads. 
Higher  fuel  loads  make  the  system  more  stable.  In  £a~*:  “-m 
higher  fuel  loads  lowered  the  stability  cutoff  inertia  ratio 

^  qv  g>  —  yw'm  q  ^  v*  ^  *s  y  ^  ^  gj  ^  **  ^  A  A  ^  fTl  ^  U  7  ^  ^  \ 

TABLE  IE.  Simulation  results  of  varying  inertia  ratio 
for  higher  fuel  loads. 


FUEL  INERTIA  RATIO 

LOAD  1.01  1.02  1.03  1.05  1.C7 

5  0%  U  M  S  S  S  S 

75%  U  M  S  S  S  S 

Note:  U  =  Unstable,  S  =  Stable,  M  =  Marginal 


The  results  of  varying  the  asymmetry  on  a  configuration  with 
an  inertia  ratio  of  1.01  and  a  fuel  load  of  75%  are  shown  in 
Table  X. 
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TABLE  X,  Simulation  results  of  varying  asymmetry  for  a 
75%  fuel  load. 


FUEL 

LOAD 

0% 

ASYMMETRY 

25% 

40% 

75% 

U 

U 

S 

Note;  U  =  Unstable,  S  =  Stable 


Here,  as  with  the  lower  fuel  loads  the  system  became  stable 
as  the  asymmetry  was  increased.  Graphing  inertia  ratio  versus 
asymmetry  for  a  configuration  with  a  75%  fuel  load  shows  that 
the  stability  cutoff  inertia  ratio  in  this  case  is  between 
1.012  and  1.015  (Figure  10).  As  with  the  symmetric  cases,  the 
stability  cutoff  inertia  ratio  is  significantly  lower  for  the 
higher  fuel  loads  than  for  the  lower  fuel  loads. 


l.U«0 

i  4 
1  02* 
1.023 
1  Z22 

1.021  4 

1.02  J 

1  Q-9 
1  c*n 


NERTIA  RATIO  VS.  %  ASYMMETRY' 

2 IS/f  M  ♦  a  1.01 


i 

A 

A 

J 

1 

J  Unstable 

U 

Stable  n  j 

j 

j 

i 

c 

□ 

O 

=  1 
i 

□ 

□ 

□ 

a 

O 

0 

,i  a 

L - r— j 

-» - ■ - 1  ■  ■  ■  . . .  ■  1  ’  - - • - 

X  A«T<VMETWY 

O  rj* 


Figure  10.  Inertia  ratio  vs.  asymmetry  for  a  75%  fuel 
1  oad . 
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3. 


ENERGY  SINK  PREDICTIONS 


As  was  shown  in  Chapter  II,  the  energy  sink  stability 
criteria  is, 

w«W"  <0  12  9) 

and , 

(Pp/Op)  +  (Pr/or)  <  0  (  -  0 ) 

where  one  or  the  other,  or  both,  of  a.  or  o„  must  be  positive 
such  that  equation  (30)  is  true.  For  the  specific  case  of  a 
dual  spin  spacecraft  with  a  d  as  pun  platform.,  aquation  (20) 

a  ^  3  ^ 


1.33,  the  I  ewer  c 

'rt"o'^ic*’icn3  w ° ** ^  '  e c^rrQr,i'  ^  o  *"  ^ 

correct  for  asymmetries  greater  t 
for  the  transition  cone  between  25' 
Spencer's  method  to  determine  the  i 
equates  to  an  inertia  ratio  of  1. 
to  an  inertia  ratio  of  1.032. 


I  — k  ->  -3  r*  V-  -o>  -v-  -»  •  *  ^  V" 

’"'rw  t  ^  -  c>c'  1  o  c  c 

n  40%,  and  indeterminate 
and  40%  asymmetry.  Using 
■rtf  a  ratio,  40%  asymmetry 
5.  25%  asymmetry  equates 
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V.  SUMMARY  AMD  CONCLUSIONS 


This  chapter  presents  the  conclusions  based  on  the  results 
presented  in  the  Chapter  IV. 

A.  SUMMARY 

The  energy  sink  stability  criteria  for  a  dual  spin 
spacecraft  with  a  despun  platform  specify  that  the  system  must 
have  an  inertia  ratio  greater  than  one.  The  results  of  this 
study  indicate  that  it  is  not  sufficient  for  the  system's 
inertia  ratio  just  to  be  greater  than  one.  It  must  be  greater 
than  one  by  a  predictable  amount.  Below  this  stability  cutoff 
inertia  ratio  the  energy  sink  method  predicts  stability  when 
the  simulation  results  indicate  the  system  is  unstable. 
Further,  the  stability  cutoff  inertia  ratio  varies  inversely 
with  the  fuel  load.  For  the  case  of  the  lower  fuel  loads  used 
in  this  study  the  stability  cutoff  inertia  ratio  is 
approximately  1.045  for  a  symmetric  spacecraft.  For  the 
higher  fuel  loads  it  is  approximately  1.02  (Figure  11). 


36 


'  i  i  1  v  *j  i  lu^'J 


n  J 

-  -I 

i  j. 
0  S3  -I 
0  9Q  -I 
'J  -J 
d  jo  -i 
■j  H5  -* 


r*jro»T-  ~  sk  ;■ 


cr  r  ^ i  i 


Inertia  ratio  vs.  fuel  load 


shewing  tor 


stability,  inertia  ratio  must  be  greater  than  one. 


As  shown  in  the  development  in  Chapter  II.  Likins’ 
stability  criteria  can  be  reduced  to  sKow  vUat  the  '<ev 
parameter  for  stability  predictions  is  the  ratio  of  the  soin 
moment  of  inertia  to  the  algebraic  mean,  transverse  moment  of 
inertia.  Recall  that  Spence'*  c1~;~ed  t**at  Likins  was 
incorrect  and  that  the  geometric  mean  transverse  moment  of 
inertia  is  the  key  stability  parameter  The  results  o*  this 
study  support  Spencer,  indicating  that  as  the  asymmetry  is 
increased  the  system  becomes  more  stable.  This  trend  is 
predictable  using  the  geometric  mean  in  computing  the  inertia 
ratio.  Interestingly,  the  stability  cutoff  inertia  ratio  for 


37 


the  asymmetric  variations  on  a  satellite  with  a  symmetric 
inertia  ratio  of  1.03  is  around  1.035.  As  the  asymmetry  is 
increased,  the  inertia  ratio  grows  from  the  symmetric  1.03 
until  it  eventually  reaches  1.05  for  a  75%  platform  asymmetry. 
It  exceeds  the  stability  cutoff  inertia  ratio  at  40% 
asymmetric  (Figure  12).  It  must  be  pointed  out  that  while 
Spencer's  geometric  mean  transverse  moment  of  inertia  more 
accurately  shows  the  effects  of  asymmetry,  it  is  also 
inaccurate  in  predicting  stability  below  the  stability  cutoff 
inertia  ratio. 


nchT iA  mATIG  Vb.  %  ASiMMETPi 


Figure  12.  Inertia  ratio  vs.  asymmetry.  Calculating 
inertia  ratio  using  Likins’  method  vs.  Spencer’s  method. 
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B.  CONCLUSIONS 


It  is  unclear  why  there  is  a  differ  ?r.c  a  *.  n  the  st  ?fc  i  1  *  *• '' 
cutoff  for  the  symmetric  3  -d  asyrret  *"  i  o  corf  igurat  i^~s  .  TUis 
is  perhaps  an  opportunity  for  future  -“search  .  ’.4 hat  has  been 
determined  or  confirmed  in  this  study  is  that: 

(1)  the  energy  sink  stability  criteria  is  not  valid  below 
a  stability  cutoff  inertia  ratio  which  is  not  Just  equal  to 
roe,  but  g-eater  than  one. 

(2)  stability  increases  as  the  fue?  1 ’.rH  •*• 


•-'■-vo  -m  o 


co-l 


/  O  N 


o*-  ;  i  i  *■ '  * 


"rooc.G'’' 


^ 0  ol  2tf  CTn 


s,ro.J 


's'  Spencer  was  correct  in  asserting  tu3t  it 
accurate  to  use  the  geometric  mean  transverse 
inertia  than  the  algebraic  mean  f-?^vD’rre  rc^er*: 

-  r o.mput i rg  the  inertia  ratio. 
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APPENDIX  A 


SIMULATION  DATA 

Included  here  are  the  graphs  of  nutation  angle  versus  time 
for  each  simulated  case.  The  graphs  ar«‘  ordered  in  the 
sequence  in  which  they  are  presented  in  the  text. 
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APPENDIX  B 


ENERGY  SINK  PREDICTIONS 

Included  in  this  appendix  are  the  spread  sheet  data  for 
the  energy  sink  predictions.  The  spread  sheet  takes  as  its 
input,  wp,  wr,  I,,  I2,  Ip,  and  Ir.  The  spread  sheet  uses  the 
following  equations  to  calculate  ar, 

h0  =  Vp  +  XrWr  (15) 

*1  =  [dp  “  I2)Wp  +  W/I,  (8) 

°2  =  [dp  -  I,)WP  +  IrWr]/I2  (9) 

a0  =  ho(Iiai  +  l2a2)/ (II  +  I2  ct2)  (20) 

ap  =  °0  "  wp  (22) 

ar  =  a0  *  (Wp  +  Wr)  (23) 

For  the  configurations  used  in  this  study,  as  long  as  ar  is 
positive,  the  energy  sink  criteria  predicts  stability. 

'"JN  '10.  'Jd  Ur  11  12  Ip  Ir  ho  Signal  Sigaa2  SigrcaO  SignaP  SigaaR  Predict  Result 

ELATION  1CTSR  (  15)  (8)  (9)  (20)  (22  )  (23) 

5TME1R1C,  Is/It  =  1.0,  753  Fill 

10030501  0.000072  3.141519  5419.4  6488.4  8012.2  6493.9  20401.29  3.138878  3.144199  3.141609  3.141536  0.000017  STABLE  191518815 
SYWETRIC,  Is/It  =  1.0,  50*.  FUEL 

90030602  0  000072  3.141519  5995.1  5984,1  7507.9  5989.6  18816.99  3.138656  3.144425  3.141611  3.141538  0.000018  STABLE  INSTABLE 
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C"N  NO.  yp  'Jr  II  12  Ip  Ir  <!tm1  Si»a2  SignaO  SionaP  Sis?aR  Predict  Result 

EQUATION  NUMBER  (IS)  (?)  (9)  (29)  (22)  '22) 

SYMMEIRIC,  Is/It  -  1.001,  75*  FUEL 

9003C79A  0.CCCO72  3.U1519  5499.4  6488.1  8018. 7  6500.4  20421  .71  3.142020  3.147347  3.144754  3.144601  0.203161  STABLE  UNSTABLE 

SYMMETRIC.  Is/It  =  1.001,  50*  FUEL 

’’2030795  0.000072  3.141519  5995.1  5984.1  7513.9  5995.6  18835.84  3.H19CO  3  147575  3.144758  3.144685  0.093165  STABLE  UNSTABLE 

SYMMETRIC,  Is/It  =  1.005,  75*  FUEL 

'329?9902  0.002972  3.141519  6499.4  6488.4  °044  7  6524.4  'O' 93  40  3.154567  3.159935  3.157332  3.157259  0.015739  STABLE  UNSTABLE 

SYMMETRIC.  Is/It  =  1.005,  50*  FUEL 

r9?2°Cl  0.09C072  3.141519  5995.1  5984.1  7537.8  6019.5  18910.92  3.154324  3.160122  3.157293  3.157221  0.915701  STABLE  UNSTABLE 

SYMMETRIC.  Is/It  =  1.01,  75*  FUEL 

92031301  0.0020/2  3.141519  6499.4  6488.4  8077.1  6550.8  20605.18  3.179249  3.175623  3.173096  3,i72933  0.031414  STABLE  UNSTABLE 

S'WETRIC,  Is/It  =  1.01,  50*  FUEL 

9CC3I292  0.000072  3.111519  '91$. i  5984.1  7567.8  6949.5  19005.17  3.170045  3.175872  3.173029  3.172956  0.031436  STABLE  UNSTABLE 

SfMMETRIC ,  Is/It  =  1.01,  26.2*  FUEL 

90021501  0.000072  3.141519  4469.5  4458.5  6026.9  4508.6  14164.29  3.169028  3.176846  3.173005  3.172932  0.031412  STABLE  UNSTABLE 

SYMMETRIC,  Is/It  =1.01,  20*  FUEL 

0.000072  3.141519  4182.5  4171.5  5737.1  4218.8  13253.86  3.168812  3.177168  3.173057  3.172984  0.031464  STABLE 

SYMMETRIC ,  Is/It  =  1.01,  15*  FUEL 

0 .000072  3.141519  3950.6  3939.6  5502.8  3984.5  12517.78  3.1685C6  3.177352  3  172995  3.172923  0.031403  STABLE 

25*  ASYMMETRIC,  Is/It  =  1.01,  75*  FUEL 

90031301  0.000072  3.141519  6872.3  6115.5  8077.1  6558.8  20605.18  2.998231  3.369256  3.173006  3.172933  0.031413  STABLE  UNSTABLE 

40*  ASYMMETRIC,  Is/It  =1.01,  75*  FUEL 

90031302  0.000072  3.141519  7096.1  5891.7  8077.1  6558.8  20605.18  2.903673  3  497237  3.173006  3.172933  0.031413  STABLE  STABLE 
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RIJN  NO.  Up  Ur  IS  !2  I?  Ir  ho  Sigsial  SigM2  SigjiaO  SignaP  SigieR  Predict  Result 
►OVATION  NUMBER  (’.51  (?)  (1)  (  20  )  (22)  (  23) 

SmiRIC,  Is/It  -1.1,  25',  F‘J5L 

0.CCCC72  3.14151?  54??. 4  6488.4  8661.6  7143.3  22441  .44  3.452776  3.458630  3.455773  3.455701  0.314181  STABLE  518815 

SYMMETRIC ,  Is/It  -1.1,  50'.  FUEL 

O.COOC72  3.141519  5995.1  5984.1  8106.9  6588.6  20698.80  3.452548  3.458894  3.455791  3.455718  0.314198  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.1,  26.2',  FUEL 

90C22001  0.000072  3.141519  1469.5  4458.5  6428.7  4910.4  15426.58  3.451451  3.159966  3.455776  3.455703  0.314183  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.1,  201  FUEL 

90022003  0.000072  3.141519  4182.5  4171.5  6113  4591.7  14134. 78  3.451161  3.460261  3.455778  3.455705  0.214185  STABLE  STABLE 

^VUWT(}TI*  T-'’f  ;  !  1  CIP 

70022002  0.000072  3.H1519  3950.6  3939.6  6857.9  4339.6  13633.36  3.45C888  3.460523  3.155771  3.455699  0.314179  STABLE  STABLE 

ft  ASYMMETRIC,  Is/It  =  1.1,  26.2',  FUEL 

90022004  0.000072  3.141519  4516.2  4411.8  6428.7  4910.4  15426.58  3.415762  3.496590  3.455776  3.455703  0.314183  STABLE  STABLE 

5».  ASYMMETRIC,  Is/It  =  1.1,  20*  FUEL 

9C0220C5  0.000072  3.141519  4227.3  4126.7  6113  4594.7  14434.78  3.414587  3.497826  3.455778  3.455705  0.314185  STABLE  STABLE 

5'.  ASYMMETRIC,  Is/It  =  1.1,  15*  FUEL 

90022101  0.000072  3.141519  3993.7  3896.5  5857.9  4339.6  13633.36  3.413647  3.498800  3.455771  3.455699  0.314179  STABLE  STABLE 

10'.  ASYMMETRIC,  Is/It  =  1.1,  26.2*  FUEL 

9C022102  0.000072  3.141519  4562.9  4365.1  6428.7  4910.4  15426.58  3.380803  3.533B98  3.455776  3.455703  0.314183  STABLE  STABLE 

10*  ASYMMETRIC,  Is/It  =1.1,  20*  FUEL 

90022103  0.000072  3.141519  4272  4082  6113  4594.7  14431.79  3.378859  3.536128  3.455778  3.4S5705  0.314185  STABLE  STABLE 

10*  ASYMMETRIC,  Is/It  =  1.1,  15*  FUEL 

90022104  0.000072  3.141519  4036.8  3853.4  5857.9  4339.6  13633.36  3.377201  3.537933  3.455771  3  455699  0.314179  STABLE  STMLE 

IS*  ASYMMETRIC,  Is/It  =  1.1,  26.2*  FUEL 

90022105  0.000072  3.141519  4609.5  4318.5  6428.7  4910.4  15426.58  3.346626  3.572131  3.455776  3.455703  0.314183  STABLE  STABLE 
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PUN  NO.  Up  Ur  11  12  !p  Ir  So  Signal  Sigma2  SigaiaO  SignaP  SignaR  Predict  Result 

EQUATION  NUMBER  (IS)  (?)  (?)  (  20)  (22)  (  23) 

IS'.  ASYMMETRIC,  Is/It  =  1.1,  20’.  FUEL 

O00221CS  0.C0C072  3.141519  1316.8  4032.2  6113  4504.7  14434.7?  3.343794  3. 575367  3.455778  3. 4S5705  0.314185  STA8LE  STABLE 

15'.  ASYMMETRIC,  Is/It  =  1.1.  151  FUEL 

90022107  0.C00072  3.141519  4080  3810.2  5857.9  4339.6  13633.36  3.341443  3.578045  3.455771  3.455699  0.314179  STABLE  STABLE 

201  ASYMMETRIC,  Is/It  =  1.1,  26.21  FUEL 

'90022108  0  000072  3.141519  4656.2  4271.8  6428.7  4910.4  15426.58  3.313061  3.611182  3.455776  3. 45S703  0.314183  STABLE  STABLE 

201  ASYMMETRIC.  Is/It  =  1.1,  201  FUEL 

90022109  0.CC9072  3.141519  4361.5  3392.5  6113  4594.7  14434.78  3.309525  3.615396  3.455778  3.455705  0.314185  STABLE  STABLE 

201  ASYMMETRIC,  Is/It  =  1.1,  151  FUEL 

90022110  O.OOOC72  3.141519  412?. 1  3767.1  5857.9  4339.6  13633. 36  3.306514  3.6189?!  3.455771  3.455699  0.314179  STABLE  ST ABLE 

251  ASYMMETRIC,  Is/It  =  1.1.  26.21  FUEL 

90022111  0.000072  3.141519  4702.9  4225.1  6428.7  4910.4  15426.58  3.280163  3.651095  3.453776  3.455703  0.314183  STABLE  STABLE 

251  ASYMMETRIC,  Is/It  =  1.1,  201  FUEL 

90022112  0.0C0072  3.141519  4406.3  3947.7  6113  4594.7  14434.78  3,275877  3.656424  3.455778  3.455705  0.314185  STABLE  STABLE 

251  ASYMMETRIC,  Is/It  =  1.1,  151  FUEL 

90022113  0.000072  3.141519  4166.2  3724  5857.9  4339.6  13633.76  3.272309  3.660865  3.455771  3.455698  0.314179  STAPLE  STABLE 

351  ASYMMETRIC,  Is/It  =  1.1,  26.21  FUEL 

90022201  0.000072  3. I4I519  4850.9  4077.1  6428.7  4910.4  15426.58  3.180088  3.783629  3.455776  3.455703  0.314183  STABLE  STABLE 

351  ASYMMETRIC,  Is/It  =  1.1,  201  FUEL 

90022202  0.000072  3.141519  4547  3807  6113  4594.7  14434.78  3,174512  3.791556  3.455777  3.455705  0.314185  STABLE  STABLE 

351  ASYMMETRIC,  Is/It  =  1.1,  151  FUEL 

90022203  0.000072  3.141519  4300.8  3589.4  5857.9  4339.6  13633.36  3.169899  3.798142  3.455771  3.455698  0.314178  STABLE  STABLE 

551  ASYMMETRIC,  Is/It  =  1.1,  26.21  FUEL 

90022204  0.000072  3.141519  4796.2  4131.8  6428.7  4910.4  15426.58  3.216355  3.733539  3.455776  3.455703  0.314183  STA8LE  STABLE 
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ptJN  “Q .  Up  'Jr  11  12  Ip  Ir  ho  Signal  Signa2  SigmaO  SigaaP  SigmaR  Predict  Result 

ECL'AMON  NUMBER  ( 15)  (8)  (?)  (20)  (22)  (  23) 

551  A5YWETRIC,  Is/It  -  1.1,  20*.  FUEL 

’0022301  0.000072  3.141519  4495.8  3858.2  6113  4594.7  14434  78  3.210664  3.741241  3.455777  3.455705  0.314185  STABLE  STABLE 

ASY>METR!C,  Is/It  -  l'l,  15%  FUEL 

’0022202  0.000072  3.141519  4252.5  3637.7  5857.9  4339.6  13633.36  3.205902  3.747713  3.455771  3.455698  0.314178  STABLE  STABLE 


SYMMETRIC ,  Is/It  =  1.03,  75%  FUEL 

90022701  0.C00072  3.141519  6499.4  6488.4  8207  6688.7  21013.28  3.233038  3.238519  3.235849  3.235776  0.094256  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.03,  501  FUEL 

90022601  0.000072  3.141519  5995.1  5984.1  7687.6  6169.3  19381.53  3.232823  3.238766  3.235865  3.235792  0.094272  STABLE  STABLE 

SYMMETRIC.  Is/It  =  1.03,  26.21  FUEL 

90022402  0.000072  3.141519  4469.5  4458.5  6116.2  4597.9  14444.83  3.231796  3.239769  3.235851  3.235778  0.094258  STABLE  UNSTABLE 

SYMMETRIC.  Is/It  =  1.03,  201  FUEL 

90022401  0.000072  3.141519  4182.5  1171.5  5820.6  4302.3  13516.18  3.231531  3.240052  3.235859  3.235786  0.094266  STABLE  UNSTABLE 

SYM74ETRIC ,  Is/It  =  1.03,  151  FUEL 

9CC22304  0.000072  3.141519  3950.6  3939.6  5581.7  4063.4  12765.65  3.231248  3.240270  3.235826  3.235753  0.C94233  STABLE  UNSTABLE 

251  ASYITCTRIC,  Is/It  =  1.03,  751  FUEL 

90022801  0.C0C072  3.141519  6856.1  6131.7  8207  6688.7  21013.28  3.064837  3.426909  3.235848  3.235776  0.094256  STABLE  STABLE 

251  ASYMMETRIC,  Is/lt  =  1.03,  SCI  FUEL 

90022802  0.000072  3.141519  6296.2  5683  7687.6  6169.3  19381.53  3.078225  3.410360  3.235864  3.235792  0.094272  STABLE  STABLE 

251  ASYMMETRIC,  Is/It  =  1.03,  26.21  FUEL 

90022501  0.000072  3.141519  4741.9  4186,1  6116.2  4597.9  14444.83  3,046149  3.450585  3.235850  3.235778  0.094258  STABLE  UNSTABLE 

251  ASYMMETRIC,  Is/It  =  1.03,  201  FUEL 

90022502  0.000072  3.141519  4442.8  3911.2  5820.6  4302.3  13516.18  3.042203  3.455681  3.235859  3.235786  0.094266  STABLE  UNSTABLE 

251  ASYM-ETRIC,  h/lt  =  1.03,  FUEL 


90022503  0.000072  3.141519  4200.7  3689.5  5581.7  4063.4  12765.65  3.038872  3.459913  3.235825  3.235753  0.094233  STABLE  UNSTABLE 
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RUN  NO.  Up  Ur  II  12 
EQUATION  NUMBER 

40*.  ASYMMETRIC,  Is/It  =  1.03,  75%  FUEL 
90022804  0.C00072  3.141519  7070.1  5917.7 

40%  ASYMMETRIC ,  Is/It  =  1.03,  50%  FUEL 
90030101  0.000072  3.141519  6476.9  5502.3 

40%  ASYMMETRIC,  Is/It  =  1.03,  26.2%  FUEL 
90030102  0.000072  3.141519  4905.4  4022.6 

40%  ASYMMETRIC,  Is/It  =  1.03,  20%  FUEL 
90020103  0.CC0072  3.141519  4599  3755 

40%  ASYMMETRIC ,  Is/It  =  1.03,  15%  FUEL 
90030104  0.000072  3.141519  4350.8  3539.4 

55%  ASYMMETRIC,  Is/It  =  1  .03,  75%  FUEL 

0.000072  3.141519  7284.1  5703.7 

55%  ASYMMETRIC,  Is/It  =  1.03,  50%  FUEL 
90022803  0.000072  3.141519  6657.6  5321.6 

55%  ASYMHETRIC,  Is/It  =  1.03,  26.2%  FUEL 
90022504  0.000072  3.141519  5068.9  3859.1 

55%  ASYMHETRIC,  Is/It  =  1.03,  20%  FUEL 
90022505  0.000072  3.141519  4755.2  3598.8 

55%  ASYMHETRIC,  Is/It  =  1.03,  15%  FUEL 
90022506  0.000072  3.141519  4500.9  3389.3 


Ip  Ir  ho  Signal  Sig«a2  SignaO  SigmaP  SignaR  Predict  Result 
(15  (8)  (9)  (20)  (22)  (23) 

8207  6688.7  21013.28  2.072072  3.550833  3.235848  3.235775  0.094255  STABLE  STABLE 

7687.6  6169.3  19381  .53  2.992347  3.522357  3.235864  3.235791  0.094272  STABLE  STABLE 

6116.2  4597.9  14444.83  2.944621  3.590832  3.235850  3.235777  0.094257  STABLE  STABLE 

5820.6  4302.3  13516.18  2.938880  3.599127  3.235858  3.235785  0.094265  STABLE  STABLE 

5581.7  4063.4  12765.65  2.934035  3.606639  3.235825  3.235752  0.094232  STABLE  STABLE 

8207  6688.7  21013.28  2.884758  3.684056  3.235818  3.235775  0.094255  STABLE 

7687.6  6169.3  19381.53  2.911131  3.641959  3.235864  3.235791  0.094271  STABLE  STABLE 

6116.2  4597.9  14444,83  2.849643  3.742963  3.235849  3.235777  0.094257  STABLE  STABLE 

5820.6  4302.3  13516.18  2.842345  3.755651  3.235857  3.235785  0.094265  STABLE  STABLE 

5581.7  4063.4  12765.65  2.836190  3.766361  3.235824  3.235752  0.094232  STABLE  STABLE 


SYMMETRIC.  Is/It  =  1.035,  26.2%  FUEL 

90030603  0.000072  3.141519  4469.5  4458.5  6138.5  4620.2  14514.89  3.247471  3.255483  3.251544  3.251472  0.109952  STABLE  UNSTABLE 

SYMMETRIC,  Is/tt  =  1.035,  20%  FUEL 

90030604  0.0C0072  3.141519  4182.5  4171.5  5841.5  4323.2  13581.84  3.247230  3.255792  3.251578  3.251505  0.109985  STABLE  UNSTABLE 
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RUN  NO.  Jp  'Jr  II  12  !p  Ir  ho  Sigmal  SigraZ  SignaO  OigoaP  SigoiaR  Predict  Result 
E3UAT1W  NUMBER  ! IS)  (8)  (0)  (  20)  (22)  (  23) 

SYMMETRIC,  Is/It  =  1.035,  15'.  Fl'El 

90030505  0.000072  3.141519  3950.6  3939.6  5601.5  4083.2  12927.86  3.246994  3.256060  3.251593  3.251520  0.110000  STABLE  UNSTABLE 

SYMMETRIC,  Is/It  =1.04,  26.21  FUEL 

'>0030701  0.000072  3.111519  4469.5  4458.5  6160.9  4642.6  [4595.26  3.263216  3.271266  3.267309  3.267236  0.125716  STABLE  UNSTABLE 

SYMMETRIC,  Is/It  =1.04,  201  FUEL 

90030702  0.000072  3.141519  4182.5  4171.5  5862.4  4344.1  13647.50  3.262928  3.271532  3.267297  3.267225  0.125705  STABLE  UNSTABLE 

SYMMETRIC,  Is/It  =  1.04,  151  FUEL 

9003C703  O.OOCC72  3.141519  3960.6  3939.6  5621.2  4102.9  12889.75  3.262659  3.271769  3. 267281  3. 2672C9  0.125688  STABLE  UNSTABLE 

SYMMETRIC,  Is/It  =  1.045,  26.21  FUEL 

'>0030303  0.000072  3.141519  1469.5  4458.5  6183.2  4664.9  14655.32  3.278890  3.286980  3.283003  3.282930  0.141410  STABLE  MARGINAL 

SYMMETRIC,  Is/It  =  1.045,  201  FUEL 

9C0308C4  0.000072  3.141519  4182.5  4171.5  5883.3  4365  13713.16  3.278627  3.287272  3.283017  3.282944  0.141424  STABLE  MARGINAL 

SYMMETRIC,  Is/It  =  1.045,  151  FUEL 

90030905  0.000072  3.141519  3950.6  3939.6  5640.9  4122.6  12951.64  3.278325  3.287479  3.282968  3.202896  0.141376  STABLE  MARGINAL 

Svr>ETRIC,  Is/It  =  1.05,  751  FUEL 

90022702  0.000072  3.141519  6499.4  6488.4  8336.9  6818.6  21421.37  3.295827  3.301415  3.298691  3.298618  0.157099  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.05,  501  FUEL 

90022602  0.000072  3.141519  5995.1  5984.1  7807.4  6289.1  19757.90  3,295602  3.301660  3.298701  3.298628  0.157108  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.05,  26.21  FUEL 

90022308  0.000072  3.141519  4469.5  4458.5  6205.5  4687.2  14725.38  3.294565  3.302693  3.298697  3.298624  0.157104  STABLE  STABLE 

SYMMETRIC ,  Is/It  =  1.05,  201  FUEL 

90022309  0.000072  3.141519  4182.5  4171.5  5904.1  4385.8  13778.50  3.294250  3.302937  3.298661  3.298588  0.157068  STABLE  STABLE 
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piN V?  Ur  II  12  Ip  !r  S  Signal  Sigr>32  Si g^aO  S i gnaP  Si g^aR  °redict  Result 

r~'(irt«t|  Vl«qr?  I  IO  (9 )  (1)  (20)  (22)  (2?) 

fT'““E,R!C,  jr/]t  t  1  09,  in*  r*tL 

’C022’C3  C.0CCC72  3.141519  3010.6  3339.6  6660.7  4142. 413013.94  3.294071  3.303263  3.298736  3. 29°(53  C. 157143  STABLE  STABLE 
253  ASYMMETRIC,  Is/It  =  1.05,  26.2%  FUEL 

’0022405  0.C00072  3.141519  4730.9  4197.2  6205.5  4687.2  14725.30  3.112597  3  508300  3.298696  3.298624  0.157104  STABLE  STAPLE 

253  ASYMMETRIC,  Is/It  =1.05,  203  FUEL 

O0C22406  0.000072  3.111519  1432.4  3021.6  ”04.1  4385.8  13778.50  3.108524  ?.51?40«  3.299660  ’.299508  0.157068  STAPLE  STAPLE 

253  ASYMMETRIC,  Ts/It  =  1.05,  153  FUEL 

’0022407  0  OCOC72  3.141519  4190.9  ’699.3  5(60.7  4142.4  12012.04  3.105103  2.51:038  3.298775  3.298663  0.157143  STABLE  STAPLE 

SYMMETRIC,  Is/It  =  1.07,  753  FUEL 

’0022703  0.CC0072  3.141519  6499.4  6188.4  8466.8  6948.5  21829.16  3.358616  3.364310  3.261534  3.361461  0.219941  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.07,  503  FUEL 

’0022603  3.000072  3.111519  59’S.l  6994.1  7927.2  6108.9  30134.26  3  358380  3.361553  3.361537  3.761464  0.219944  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.07,  26.23  FUEL 

’0022103  0.0C0072  3.141519  4469.5  4458.5  6294.8  4776.5  15005.92  3.357333  3.365616  3.361542  3.361470  0.219950  STABLE  STA8LE 

SYMMETRIC,  Is/It  =  1.07,  203  FUEL 

’0022104  0.000072  3.111519  4182.5  4171.5  5907.7  4469.4  14041.14  3.357045  3.365897  3.761538  3.961465  0.219945  S’'°LE  STABLE 

SYMMETRIC,  Is/It  =  1.07,  153  FUEL 

’0022705  0.000072  3.141519  3950.6  3939.6  5739.6  4221.7  13261.71  3.356813  3.366186  3.361566  3.361493  0.219973  STAPLE  STABLE 

<Y>~ETR!C,  Is/It  =  1.02,  753  FUEL 

’0073703  0.000072  3.111519  6499.4  6488.4  9142.1  6623.8  20809,39  3.201667  3.207095  3.204452  3.204379  0.062859  STABLE  STABLE 

SYMMETRIC,  Is/It  =  1.02,  503  FUEL 

90022*03  0.000072  3.141519  5995.1  5984.1  7627.7  6109.4  19193.35  3.201434  3.207319  3.204447  3.204374  0.062854  STAPLE  STABLE 
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